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Our knowledge of the mass composition of cosmic rays is deficient at all energies above 10 17 eV. Here, systematic 
differences between different measurements are discussed and, in particular, it is argued that there is no compelling 
evidence to support the common assumption that vast majority of the cosmic rays of the highest energies are 
protons. Our knowledge of the mass needs to be improved if we are to resolve uncertainties about the energy 
spectrum. Improvement is also needed for proper interpretation of data on the arrival direction distribution of 
cosmic ray. 



1. The Scientific Motivation for Studying 
the Highest Energy Cosmic Rays 

Since the recognition in 1966, by Greisen and 
by Zatsepin and Kuzmin, that protons with en- 
ergies above 4 x 10 19 eV would interact with the 
cosmic microwave radiation, there has been great 
interest in measuring the spectrum, arrival di- 
rection distribution and mass composition of ul- 
tra high-energy cosmic rays (UHECR), defined as 
those cosmic rays having energies above 10 19 eV. 
Specifically they pointed out that if the sources 
of the highest energy protons are universally dis- 
tributed, there should be a steepening of the en- 
ergy spectrum in the range from 4 to 10 x 10 19 eV. 
This feature has become known as the GZK 'cut- 
off' but the sharpness of the steepening expected 
depends on unknown factors such as the evolu- 
tion and production spectrum of the sources. If 
the UHECR are mainly Fe nuclei then the spec- 
trum is also expected to steepen, but it is harder 
to predict the character of this feature as the rele- 
vant diffuse infrared photon field is poorly known: 
steepening is expected to set in at higher energy. 

2. Importance of mass composition for ac- 
curate estimates of the energy spectrum 

The energy spectrum of cosmic rays has fre- 
quently been inferred from observations of sig- 
nals in arrays of scintillators or water-Cherenkov 
detectors (see Nagano and Watson 1, for a re- 



view). In practice, the energy is derived from a 
measurement of the detector response at (typi- 
cally) 600 m from the shower axis using the re- 
sults of detailed Monte Carlo calculations. For ex- 
ample, in the most recent reports of spectra from 
the AG AS A [2] and Haverah Park groups [3], 
the QGSJET model of high energy interactions 
has been adopted. Additionally, in the AGASA 
work, a study of the impact of different models 
was made and it was found, for example, using 
QGS JET98 in the Corsika propagation code, that 
the energy estimates for protons or iron nuclei 
differ by about 12% near 10 20 eV (with protons 
giving the higher energy), while with SIBYLL 1.6 
the difference is 19%, in the same direction and 
the overall energy estimates are about 5% higher. 
These two sources of systematic error in energy 
estimates — from mass and model — are essential 
to keep in mind when comparing data sets. As 
the relevant centre of mass energy at 10 20 eV is 
well beyond that anticipated at the LHC, it is 
clear that the systematic effect in the model is 
extremely difficult to quantify. 

The other method that has been used to mea- 
sure primary energies is the fluorescence tech- 
nique in which photons emitted by excited N2 are 
detected. With this method, it is possible to make 
a calorimetric estimate of the large fraction of the 
primary energy that is transferred into ionisation 
as the shower particles cross the atmosphere. The 
estimate of this part of the primary energy loss 
is relatively model independent with, as first dis- 
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cussed by Linsley @], and more recently by Song 
et al. 0, a correction of around 10% being re- 
quired for energy that does not go into atmo- 
spheric ionisation. The most recent journal pre- 
sentation of a fluorescence spectrum by the HiRcs 
group is one in which monocular data from two 
detectors are reported Note that the correc- 
tion for missing energy is mass dependent, with 
an energy estimate about 5% higher being made 
if the primaries are assumed to be Fe nuclei. 

When making comparisons of spectra, a log-log 
plot of J (the differential intensity) against en- 
ergy, E, has significant advantages over the JE 3 
vs. E plots commonly adopted in recent years. 
In particular, propagators of the latter style (in- 
cluding the present author) almost always ignore 
the fact that the error bars on such a plot should 
be shown as diagonal lines and that the uncer- 
tainty in energy (often justifiably omitted in the 
x-direction because of the bin size) should be 
added in quadrature to the uncertainty in inten- 
sity, as it comes in as the third power. The JE 3 
vs. E representations do the data a disservice 
as the incorrect error assignments serve to give 
an erroneous impression of the level of incompat- 
ibility between the data sets. For example, the 
differences between the AGASA and HiRes spec- 
tra could largely be resolved, except perhaps at 
the very highest energies, if one or other of the 
energy estimates had a systematic error of only 
~ 20 — 30%. Some movement towards reconcilia- 
tion would arise if the primary particles were iron 
nuclei. 



3. The mass of UHECR 

Knowledge about the mass of primary cosmic 
rays at energies above 10 17 eV is rudimentary. 
Different methods of measuring the mass give dif- 
ferent answers and the conclusions are dependent, 
in all cases, upon the model calculations that are 
assumed. Results from some of the techniques 
that have been used are now described and the 
conclusions drawn reviewed. Some of these tech- 
niques will be applicable with the Pierre Auger 
Observatory [7|. 



3.1. The Elongation Rate 

The elongation rate is the term used to describe 
the rate of change of depth of shower maximum 
with primary energy. The term was introduced 
by Linsley 8 and, although his original approach 
has been, to some extent, superseded by the re- 
sults of detailed Monte Carlo studies, the concept 
remains useful for organising and describing data 
and calculations. In figurenthere is a summary of 
measurements of the depth of maximum together 
with predictions from a variety of model calcula- 
tions reproduced from It is clear that, if cer- 
tain models are correct, one might infer that the 
primaries above 10 19 eV are dominantly protons 
but that others are indicative of a mixed compo- 
sition. In particular, the QGSJET set of models 
(basic QGSJET01 and the 5 options discussed 
in [H]) and the Sibyll 2.1 model force contrary 
conclusions. Note that the experimental data on 
Xmaz do not yet extend beyond 25 EeV. 
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Figure 1. The depth of maximum, as predicted 
using various models, compared with measure- 
ments. The predictions of the five modifications 
of QGSJET, discussed in [9], from which this di- 
agram is taken, lie below the dashed line that 
indicates the predictions of QGSJET01. 
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3.2. Fluctuations in Depth of Maximum 

Further insight is expected to come from the 
magnitude of the fluctuation of the position of 
depth of maximum, X maa; . If a group of showers 
is selected within a narrow range of energies, then 
fluctuations about the mean of X mol (or in a pa- 
rameter that is closely related to X moa; such as 
the steepness of the lateral distribution function 
or the spread of muon densities) are expected to 
be larger for protons than for iron nuclei. A re- 
cent study applying this idea to X maa; has been 
reported by the HiRes group ^Oj for 553 events 
above 10 EeV. 

In it is argued that the fluctuations are 
so large that a large fraction of light nuclei must 
be present in the primary beam (figure 0) : this 
conclusion is independent of whether Sibyll or 
QGSJET models are used. Clearly an under- 
standing of the tails on the high-X ma:c side of 
the histograms is crucial to the conclusion that 
they are best described by 'a predominantly light 
composition'. Of the 553 events in the data 
set, detailed atmospheric information is available 
for 419: the rest (134) occurred during 'good 
weather' and were analysed assuming average at- 
mospheric conditions. This assumption may not 
represent reality, as it is possible that the at- 
mosphere deviates from the standard conditions 
from night to night and even during a night of 
observation. This view is strengthened by the re- 
sults of balloon flights made from Malargue |11) . 
which have shown that the atmosphere changes in 
a significant way both diurnally and seasonally. If 
a standard atmosphere is used, some of the fluc- 
tuations observed in X„ la:E may be attributed in- 
correctly to shower, rather than to atmospheric, 
variations. It would be instructive to see the dis- 
tributions in X maa: for the two weather groups. 

Questions about the interpretation are also 
raised when the Monte Carlo analysis of uncer- 
tainties in X maa; and the fluctuations in X max are 
considered (Perrone ^2]). He finds that at dis- 
tances beyond 20 km, there are significant sys- 
tematic shifts in the X maa; values derived and in 
the spread of the X max values. At 20 km, X maa; 
is estimated to be 60 g cm -2 deeper in the atmo- 
sphere, on average, than reality and the fluctua- 
tions in X„ ia:r for iron nuclei are considerable with 



a ~ 100 g cm~ 2 . These factors act in such a way 
as to suggest that the elongation rate reported 
by HiRes (and presumably also by Fly's Eye) may 
have been systematically over-estimated and that 
the fluctuations in X maa; are not due entirely to 
protons. Thus, it may be premature to draw con- 
clusions about the presence of light nuclei from 
the analyses of fluctuations as presented so far. 

A further issue of some concern is the quality 
of events selected for the HiRes analysis. In ^U] 
it is stated that events were selected when x 2 per 
degree of freedom for the fit of the longitudinal 
development was less than 20. This is a rather 
loose cut. Furthermore, the uncertainties in X mQa; 
are required to be less than 200 g cm -2 for the fit 
made using both Eyes (as compared with 400 g 
cm" 2 for each Eye). It is hard to reconcile these 
numbers with the resolution of 30 g cm -2 claimed 
for the measurements from Monte Carlo studies 
and it seems reasonable to question whether all 
of the events in the tails in the data of figure [2] 
arise from the presence of light primary particles. 

3.3. Mass from muon density measure- 
ments 

It is well known that a shower produced by 
an iron nucleus will contain a greater fraction 
of muons at the observation level than a shower 
of the same energy created by a proton primary. 
Many efforts to derive the mass spectrum of cos- 
mic rays have been attempted using this fact over 
the full range of air shower observations. How- 
ever, although the differences are predicted to 
be relatively large (on average there are ^70% 
more muons in an iron event than a proton event), 
there are large fluctuations and, again, there are 
differences between what is predicted by partic- 
ular models. Thus, the QGSJET model set pre- 
dicts more muons than the Sibyll family, the dif- 
ference arising from different predictions as to the 
pion multiplicities produced in nucleon-nucleus 
and pion- nucleus collisions that in turn arise from 
differences in the assumptions about the par- 
ton distribution within the nucleon |13| . In a 
contribution to these proceedings, Shinosaki has 
described the data on muons signals from the 
AGASA array. There are 129 events above 10 19 
eV, of which 19 have energies greater than 3x 10 19 
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Figure 2. Presentation in [10] of the distribution of the measured values of X max for the energy intervals 
1 to 2.5 EeV and 2.5 to 25 EeV. The higher energy range is in the lower plots and the data (solid lines) 
are compared with the QGSJET model (dashed) and the Sibyll model (dotted). 



eV. Measurements of muon densities at distances 
between 800 and 1600 m were used to derive the 
muon density at 1000 m with an average accu- 
racy of 40%. This quantity is compared with the 
predictions of model calculations. The difference 
between the proton and iron predictions is small, 
especially when fluctuations are considered. The 
AGASA group conclude that above 10 19 eV the 
fraction of Fe nuclei is < 40% at the 90% confi- 
dence level. In my view, the 5 events above 10 20 
eV for which such measurements are possible, are 
fitted as well by iron nuclei as by protons. Fur- 
ther, the conclusions are sensitive to the model 
used: as the Sibyll model predicts fewer muons 
than the QGSJET model, higher iron fractions 
would have been inferred had that model been 
adopted. At this meeting, Ostapchenko has dis- 
cussed changes to the QGSJET model that will 
reduce the number of muons expected for a pro- 
ton primary (see also Engel |14|V The magnitude 
of the effect is not yet clear but it is in such a di- 
rection as to raise the fraction of iron nuclei at 
the energies in question. 

At lower energies, there are muon data from 
the Akeno array and from AGASA [T5\ . Different 
analyses have been made of these. The AGASA 



group claim that the measurements are consistent 
with a mass composition that is unchanging be- 
tween 10 18 and 10 19 eV. Another interpretation 
is discussed below. 

3.4. Mass estimates from the lateral distri- 
bution function 

The rate of fall of particle density with distance 
from the shower axis provides a further parameter 
that can be used to deduce the mass composition. 
Showers with steeper lateral distribution func- 
tions (LDFs) than average will arise from showers 
that develop later in the atmosphere, and vice 
versa. A detailed measurement of the LDFs of 
showers produced by primaries of energy greater 
than 10 17 eV was made at Haverah Park using a 
specially constructed 'infilled array' in which 30 
additional water tanks of 1 m 2 were added on a 
grid with spacing of 150 m. When the experimen- 
tal work was completed in 1978, the data could 
not be fitted with the shower models then avail- 
able for any reasonable assumption about the pri- 
mary mass. Recently JH], these data have been 
re-examined using the QGSJET98 model. The 
choice of model was justified by showing that it 
adequately described data on the time spread of 
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the Haverah Park detector signal over a range of 
zenith angles and distances near the core (< 500 
m). Here the difference predicted between the 
average proton and iron shower is only a few 
nanoseconds and the fit achieved is good. Density 
data were fitted by a function p(r) - r (v+r/iooa) ^ 
where r\ is the steepness parameter. The spread 
of rj was compared with predictions for different 
primary masses. The proton fraction, assuming a 
proton-iron mixture, is found to be independent 
of energy in the range 3 x 10 17 to 10 18 eV and 
is (34 ± 2) %. When this fraction is evaluated 
with QGSJET01, in which a different treatment 
of diffractive processes is adopted from that in 
QGSJET98, then the fraction rises to 48%. The 
fraction is larger because the later model predicts 
shower maxima that are higher in the atmosphere 
and accordingly, to match the observed fluctua- 
tions, the proton fraction must be increased. The 
difference in the deduced ratio thus has a system- 
atic uncertainty from the models that is larger 
than the statistical uncertainty. Although the 
necessary analysis has not been made, it is clear 
that the Sibyll 2.1 model would be consistent with 
an even smaller fraction of protons. 

A similar analysis has been carried out using 
data from the Volcano Ranch array. As with the 
Haverah Park information, no satisfactory inter- 
pretative analysis was possible when the mea- 
surements were completed. Using 366 events 
for which Linslcy left detailed information, and 
QGSJET98, the fraction of protons is estimated 
as (11 ± 5(stat) ± 12 (sys))% at ~ 1 EeV [12]. 
With QGSJET01 the flux of protons would in- 
crease to ~ 25%, indicating again that model un- 
certainties remain a serious barrier to interpreta- 
tion. 

3.5. Mass from the thickness of the shower 
disk 

The particles in the shower disc do not arrive 
at a detector simultaneously, even on the shower 
axis. The arrival times are spread out because of 
geometrical effects, velocity differences, and be- 
cause of delays caused by multiple scattering and 
geomagnetic deflections. The first particles to ar- 
rive (except very close to the shower axis) are the 
muons as they are scattered rather little and geo- 



metrical effects dominate. At Haverah Park four 
detectors, each of 34 m 2 , provided a useful tool 
for studying how the thickness of the shower disc 
depends upon the development of the cascade. 
Recently, an analysis of 100 events of mean en- 
ergy ~ 10 19 eV has shown that the magnitude 
of the risetime is indicative of a large fraction 
(~80%) iron nuclei at this energy ^B]- This type 
of study will be considerably extended with the 
Pierre Auger Observatory in which the photomul- 
tipliers within each water tank are equipped with 
25 ns flash ADCs. 

3.6. Summary of data on primary mass 
above 0.3 EeV 

In figure taken from , the results taken 
from various reports of the Fe fraction are shown. 
It is disappointing that the data from Volcano 
Ranch and from Haverah Park are not in better 
agreement as a similar quantity, the lateral dis- 
tribution function of the showers, was measured 
at each array and the same model - QGSJET98 
- was used for interpretation, although with dif- 
ferent propagation codes (AIRES and CORSIKA 
respectively). We cannot explain this difference: 
at 10 18 eV the estimates of the fraction of Fe are 
separated by over 2 standard deviations. 

In figure |3 data from the Akeno/AGASA and 
the Fly's Eye experiments are also shown. The 
Akeno/AGASA groups measured the muon den- 
sities in showers, normalised at 600 m. The en- 
ergy thresholds for Akeno and AGASA were 1 
and 0.5 GeV respectively. The Fly's Eye data 
are deduced from measurements of the depth of 
shower maximum. In an effort to reconcile dif- 
fering claims made by the two groups of the 
trend of mass composition with energy, Dawson 
et al. J2| reassessed the situation using a single 
model, SIBYLL 1.5 on both data sets. SIBYLL 
1.5 was an early version of the SIBYLL family 
that evolved to SIBYLL 1.6 and 1.7. It is the 
estimates of the Fe fractions from that are 
shown in figure [3] There are major discrepancies 
between these estimates and between those from 
Volcano Ranch and Haverah Park. However, the 
predictions of the muon density and of the depth 
of shower maximum made with the version of 
SIBYLL used in J5j differ significantly from those 
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Figure 3. The fraction of Fe nuclei as a function 
of energy as reported from various experiments. 
Fly's Eye (A), AGASA A100 (■), AGASA Al 
(□), Haverah Park (O) and Volcano Ranch (•) 
(Figure taken from [17]). 



that would be derived now using QGSJET98 or 
01 (or with the later SIBYLL version, 2.1). We 
now discuss this point in some detail. 

An extremely useful set of comparisons of the 
predictions from SIBYLL 1.7 and 2.1 with those 
from QGSJET98 has been given in 13 . We un- 
derstand that SIBYLL 1.6 and SIBYLL 1.7 differ 
only in that the neutral pions were allowed to in- 
teract in the latter model and it is not believed 
that this will make a serious difference to the pre- 
dictions at energies below 10 19 eV (23). There- 
fore, in what follows, we regard the SIBYLL 1.7 
and the QGSJET98 differences as being identi- 
cal to those that exist between SIBYLL 1.6 (or 
1.5) and QGSJET98, for which no similar com- 
parisons are available. It is convenient to com- 
pare conclusions at 10 18 eV. More detailed cross- 
checks, over a range of energies, would require 
more extensive knowledge of features of the Fly's 
Eye and Akeno/ AGASA systems than we possess. 



Turning first to the data from the depth of 
maximum, we note that at 10 18 eV the measured 
value of X ma x is ~675 g cm~ 2 , with an error 
that is less than the size of the data point (< 
10 g cm -2 ). The predictions for proton primaries 
made with SIBYLL 1.7 and QGSJET98 are 760 
and 730 g cm -2 respectively jlSj. Thus, a mass 
composition less dominated by Fe is favoured 
compared with the ~ 90% estimated in P3| . The 
choice of SIBYLL 2.1 would alter this argument 
rather little as the predicted depth at 10 18 eV is 
740 g cm -2 |13| . Further study of this matter 
could be made but the data from Fly's Eye has 
now been complemented by data from the HiRes 
stereo system ^U] (although the differences be- 
tween the data sets have not been explained) and 
there will also be data from the Auger instru- 
ments. 

A qualitative statement about the shift ex- 
pected in the Fe fraction, as estimated from the 
measurement of muon densities at 600 m, from 
changes in the model can be made using infor- 
mation in JSj- Although the calculations do not 
exactly match the energies of the Akeno/ AG AS A 
measurements (> 0.3 GeV is computed and > 
0.5 GeV measured), ratios between the predic- 
tions of different models are not strongly depen- 
dent upon energy threshold. What is of impor- 
tance is the ratio of the number of muons pre- 
dicted, at 10 18 eV, for SIBYLL 1.7, SIBYLL 2.1 
and QGSJET98. At 10 18 eV, the numbers arc in 
the ratios 1: 1.17: 1.44. The difference in muon 
number between SIBYLL 1.7 and QGSJET98 is 
comparable to that expected between proton and 
Fe primaries (~50%, but also model dependent). 
It is clear that the more recent models, if applied 
to the Akeno/ AG AS A data after the manner of 
the analysis of [T§|, would lead to a significant 
reduction in the predicted fraction of Fe nuclei. 
To pursue this further would require knowledge 
of the predicted densities at 600 m, information 
that is presently lacking. We note that the shift 
in the Fe fraction from the muon data is proba- 
bly substantially larger than it is when using the 
data on X max . 

We are not able to use the information re- 
ported from the HiRes-MIA experiment [2] in 
which muons and X maa; were observed simultane- 
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ously. As with Akeno/AGASA, the muon density 
at 600 m was determined. The problem we have 
is that while the papers describe the data as being 
consistent with a mass composition that becomes 
lighter with energy, this appears, on close scrutiny 
of figures ^ and [21 of |21|. to be true only for the 
X maa; data. The muon data, which are compared 
with predictions of QGSJET98, look to be consis- 
tent with a constant and heavy mass from 5 x 10 16 
to beyond 10 18 eV. It would be very interesting 
to establish that the same model gives different 
predictions for the mass variation with energy for 
different measured quantities: this might lead to 
further understanding of the models, or of the 
systematic errors in measurements of X ma;r , as 
discussed above. 

This discussion is intended to demonstrate the 
difficulties with which one is faced with when try- 
ing to compare data. Measurements from differ- 
ent experiments are rarely analysed contempora- 
neously and the shifts in the inferences from the 
use of different models can be substantial. 

4. Possibilities of Identifying Photon and 
Neutrino Primaries 

4.1. Super- heavy relic particles 

An idea to explain the UHECRs that have been 
reported beyond 100 EeV is that super-heavy 
relic particles with masses of ~ 10 12 GeV, pro- 
duced in the early Universe, may decay to pro- 
duce UHECR 122- While details of the fragmen- 
tation of these particles remains a matter of de- 
bate, it is generally accepted that the resulting 
UHECR beam would contain copious fluxes of 
neutrinos and photons. 

4.2. Limits to the fraction of photon pri- 
maries 

It is unlikely that the majority of the events 
claimed to be near 10 20 eV have photons as par- 
ents as some of the showers seem to have nor- 
mal numbers of muons, the tracers of primaries 
that are nuclei a , (see paper in these proceedings 
by Shinosaki). It has been argued that the cas- 
cade profile of the most energetic Fly's Eye event 

"This assumes that the photo-pion production cross- 
section behaves 'normally' 



|2"3] is inconsistent with that of a photon pri- 
mary |24| . However, in a recent paper by Risse 
et al. |2S] (and Risse, these proceedings) have 
used calculations made with the QGSJET01 and 
SIBYLL2.1 models to show that the profile of this 
event can be explained under the assumption of 
any baryonic primary between a proton and iron 
nucleus, and that the primary photon hypothesis, 
although not favoured, cannot be rejected. 

An alternative method of searching for pho- 
tons has been developed using showers incident 
at very large zenith angles. Deep-water tanks 
have a good response to such events out to be- 
yond 80°. At such angles the bulk of the showers 
detected are created by baryonic primaries but 
they are distinctive in that the electromagnetic 
cascade stemming from neutral pions has been 
almost completely suppressed by the extra thick- 
ness of atmosphere penetrated. At 80° the at- 
mospheric thickness is ~ 5.7 atmospheres at sea- 
level. At Haverah Park, showers at such large 
zenith angles were observed and the shower disc 
was found to have a very small time spread. A 
complication for the study of inclined showers is 
that the muons, in their long traversal of the at- 
mosphere, are very significantly bent by the ge- 
omagnetic field. A study of this has been made 
and it has been shown that the rate of trigger- 
ing of the Haverah Park array at large angles can 
be predicted ln addition, it was found that 
the energy of the primaries could be estimated 
with reasonable precision so that an energy spec- 
trum could be derived. The concept of using the 
known, and mass independent, spectrum deduced 
from fluorescence detectors to predict the trigger- 
ing rate as a function of the mass of the primary 
has led to a demonstration that the photon flux 
at 10 19 eV is less than 40% of the baryonic com- 
ponent a conclusion similar to that of the 
AG AS A group, made by searching for showers 
which have significantly fewer muons than nor- 
mal EE]- 

The flux of photons expected at 10 19 eV from 
super-heavy relic particles has recently been re- 
assessed and is expected to be lower than 
originally predicted. This arises because of a 
more detailed examination of the fragmentation 
process. A large flux of photons is now predicted 
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at 100 EeV and can be sought in the Auger data 
using the method developed at Haverah Park. 

4.3. Neutrino Primaries: 

Neutrino primaries may be detectable by 
studying very inclined showers. This idea was 
first proposed by Berezinsky and Smirnov |30| 
and was re-examined in the context of the Auger 
Observatory by Capelle et al. [3T]. A discussion 
of the potential of the Pierre Auger Observatory 
to detect such events is given in the paper by K- 
H Kampert in the proceedings of this meeting. A 
neutrino can interact anywhere in the atmosphere 
with equal probability. 

5. The Zatsepin effect 

One method capable of giving an estimate of 
the primary mass, that makes no assumptions 
about the particle physics at extreme energies, 
is that proposed by Zatsepin [121 m 1951. He 
pointed out that heavy nuclei would undergo pho- 
todisintegration in the radiation field of the sun. 
This is the same process, but with the 2.7 K and 
IR radiation fields, that would affect the spectrum 
of Fe nuclei at the very highest energies. The 
process in the solar photon field is important for 
56 Fe nuclei at energies of around 10 18 eV and the 
resulting fragments ( 55 Fe and a neutron) would 
produce showers at roughly the same time and 
at nearly the same place in the atmosphere. In 
a later calculation, Gerasimova and Zatsepin |32| 
overlooked the effect of the interplanetary mag- 
netic field and predicted a separation of fragments 
of only a few 100 metres. When the effect of the 
interplanetary field was included the core separa- 
tion was estimated to be many kilometres' 3 . The 
importance of Zatsepin's idea is that the ratio of 
the sizes of the two correlated showers would be 
proportional to the mass of the fragmented pri- 
mary. No assumptions about particle physics are 
needed. 

In view of the scale of the Pierre Auger Ob- 
servatory (~ 31 km radius) and the improved 
knowledge of the interplanetary field, a further 

^During the Pylos meeting, G T Zatsepin credited M 
Shapiro with pointing out the effect of the interplanetary 
magnetic field in discussions during the 5th International 
Cosmic Ray Conference (Guanajauto 1955) 



study was made of this effect [33]. However, it 
was found that core separations of less than 10 
km (a scale relevant to the Auger surface array 
and to AGASA) were infrequent (0.3 per year on 
Auger) and hard to detect, as the optimum en- 
ergy of about 6 x 10 17 eV is rather low. However, 
there is a stronger signal for separations of ~ 10 3 
km, the approximate distance between two sites 
being considered for the northern Auger Obser- 
vatory in Utah and South Eastern Colorado: this 
prospect may be worthy of more detailed scrutiny. 
A model-independent measurement of the mass 
of primary cosmic rays is highly desirable. This 
beautiful idea may provide a solution. 

Conclusions 

To make full use of forthcoming information 
on the energy spectrum and arrival direction dis- 
tribution at the highest energies, and to interpret 
what already exists, it is necessary to improve our 
knowledge of the mass of the cosmic rays above 
10 19 eV. Such evidence as there is does not sup- 
port the widely adopted assumption that all of 
these cosmic rays are protons: there may be a 
substantial fraction of iron nuclei present, even 
at 10 20 eV. Photons do not appear to dominate 
at 10 19 eV. 
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